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Abstract Where spatial gradients in the ampUtude of an Alfven wave are non- 
zero, a nonlinear magnetic-pressure gradient acts upon the medium (commonly 
referred to as the ponderomotive force). We investigate the nature of such a force 
in inhomogeneous 2.5D MHD plasmas by analysing source terms in the nonlinear 
wave equations for the general case of inhomogeneous B and p, and consider 
supporting nonlinear numerical simulations. Our equations indicate there are 
two distinct classes of ponderomotive effect induced by Alfven waves in general 
2.5D MHD, each with both a longitudinal and transverse manifestation; 

i) Geometric Effects: Gradients in the pulse geometry relative to the background 
magnetic field cause the wave to sustain cospatial disturbances, the longitudinal 
and transverse daughter disturbances - where we report on the transverse dis- 
turbance for the first time. 

ii) V(c^) Effects: Where a pulse propagates through an inhomogeneous region 
(where the non-zero gradients in the Alfven-speed profile ca are non-zero), the 
nonlinear magnetic-pressure gradient acts to accelerate the plasma. Transverse 
gradients (phase mixing regions) excite independently propagating fast magne- 
toacoustic waves (generalising the result of Nakariakov et al. 1997, Solar Physics, 
175, 93) and longitudinal gradients (longitudinally dispersive regions) perturb 
along the field (thus creating static disturbances in /3 = 0, and slow waves in 

We additionally demonstrate that mode conversion due the nonlinear Lorentz 
force is a one-way process, and does not act as a mechanism to nonlinearly 
generate Alfven waves due to propagating magnetoacoustic waves. We conclude 
that these ponderomotive effects are induced by an Alfven wave propagating in 
any MHD medium, and have the potential to have significant consequences on 
the dynamics of energy transport and aspects of dissipation provided the system 
is sufficiently nonlinear and inhomogeneous. 
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1. Introduction 

Due to the abundance of observational data confirming the existence and ubiq- 
uity of MHD wave motions in coronal plasma, it is clear that a well-developed 
theory of MHD wave propagation is required for understanding many dynamic 
processes ongoing in the coronal plasma (see reviews by, e.g., De Moortel 
Nakariakov & Verwichte ,2005; R uderman & Erdelyi 2009 Goossens et al. 
De Moortel & Nakariakov 2012). The magnetic topology of the corona supports 



2005 



2011 



a wide variety of features, and as such can be highly inhomogeneous. Since the 
behaviour of MHD waves is determined by the medium in which they propagate, 
an understanding of wave propagation in inhomogeneous media is essential. 

One important consequence of considering propagation in an inhomogeneous 
medium is the concept of phase mixing. Simply, as a pulse or wave train propa- 
gates through the inhomogeneous medium, oscillations on neighbouring fieldlines 
become out of phase as local Alfven speeds on each fieldline differ. This process 
naturally creates transverse gradients in the wave across the background mag- 
netic field, and acts as a (linear) dissipative mechanism of Alfven wave energy, 
and hence is often considered as a potential mechanism of coronal heating (see, 
e.g., Heyvaerts & Priest [T9831 B rownmg [T99T| Narain & Ulmschneider [T990l 



1996). Typically, such studies have considered an inhomogeneity in density (to 



permit a unidirectional field, which makes mathematics more tractable), however 
as this paper will re-enforce, the process of phase mixing is strictly dependent on 
a non-uniform background Alfven speed, i.e. is dependent on both density and 
magnetic field profiles, and is not restricted to cases of density inhomogeneity 
only. 

There is a secondary consequence of these non-zero spatial gradients in the 
amplitude of a propagating MHD wave: nonlinear magnetic-pressure gradients 
and nonlinear magnetic tension (i.e. the nonlinear Lorentz force) act upon the 
plasma. In the MHD context, these forces are often referred to as the ponderomo- 
tive force. Generally, a ponderomotive force is defined as a basic nonlinear force 
consisting of spatial gradients in a wave-field which has a non-vanishing effect 
when averaged over the period, however we note that the term has been used to 



refer to other forces (see Allan et al. 1991 and Verwichte 1999 for discussions of 
the terms historical usage) . The concept exists in areas of plasma physics other 
than MHD waves, in particular in the study of laser-plasma interaction where 



the force causes, including self-focusing (e.g. Chen 1984) and channel formation 



e.g. Boyd & Sanderson 2003). The ponderomotive force in the MHD waves 



context (i.e. the nonlinear Lorentz force arising in fluid equations of motion) has 



been formally discussed by Dewar( 1970 1 and Webb et al. ( 2005 1 . For a thorough 
discussion of the ponderomotive force in various contexts, , and a demonstration 
of the equivalence of particle-orbit ponderomotive force, see Allan et al. ( 1991 
section 3 and appendix A). 

In our case (propagating MHD waves in inhomogeneous media), the effects 
of such a force are primarily of interest in that they have the potential to 
facilitate nonlinear mode conversion, the result of which has consequences on 
energy transport and dissipation. In many previous studies of ponderomotive 
effects (e.g., in almost all of the forthcoming references), because such studies 
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are concerned with phenomena due to propagating Alfven waves, the authors 
use the term ponderomotive force to refer specifically to the nonlinear magnetic 
pressure force (as Alfven waves do not act on the medium via nonlinear magnetic 
tension, e.g. as shown in fjs]). However, as reported in Thurgood & McLaughlin 
( 2012 ) , propagating magnetoacoustic waves have a similar action on the medium, 



but in their case due to both nonlinear magnetic tension and nonlinear magnetic 
pressure. To avoid ambiguity, henceforth when referring to the force responsible 
for the investigated phenomena we avoid the term 'ponderomotive force' in 
favour of the more specific 'nonlinear magnetic-pressure gradient' or 'nonlinear 
magnetic tension' as appropriate. We use the term 'ponderomotive effect' to 
refer to the wider family of phenomena caused by the nonlinear Lorentz force 
of propagating MHD waves, and in this paper we primarily consider the specific 
set of ponderomotive effects induced by propagating Alfven waves. Outside of 
the coronal plasma context, ponderomotive effects of MHD waves have been 



investigated in magnetospheric physics (e.g. Rankin et al. 1994 
Tikhonchuck et al. 1995; Allan & Manuel 1996) and as an acceleration mech- 
anism in the solar wind (eg. Stark et al. 1995). Within the context of coronal 



plasmas, ponderomotive effects where considered initially by Hollweg (1971) and 
later by Nakariakov et al. ( |1997[ ), Verwichte (1999) and Verwichte et al. ( |1999[ ), 
the last three of which are key references in this paper. 



Verwichte ( 1999 1 and Verwichte et al. ( 1999 ) considered the nonlinear evo- 



lution of an Alfven wave in a homogeneous 1.5D model. They found that the 
nonlinear magnetic-pressure gradient was responsible for two features of the 
resultant wave behaviour. Firstly, they found that the propagating Alfven wave 
sustains a disturbance which is longitudinal to the background magnetic field, 
and caused no net perturbation to the plasma (since the average nonlinear 
magnetic pressure over the wave period was zero). They dubbed this feature 
a ponderomotive wing. Secondly, they found that the interaction between two 
crossing Alfven waves gives rise to a non-zero net perturbation of plasma den- 
sity. Hence, in warm plasmas, this cross-ponderomotive effect was identified 
as a potential mode conversion mechanism, as would certain combinations of 
pulse geometry and background inhomogeneity leading to non-zero longitudinal 
perturbations as the pulse propagates. 



Nakariakov et al. (1997) considered the nonlinear excitation of the fast wave 



by a propagating Alfven wave which undergoes phase mixing in a unidirectional 
field structured by a background density profile (hence, an inhomogeneous back- 
ground Alfven speed). They derive governing wave equations, then evaluate them 
at an initial instant where only a pure, linear Alfven wave is perturbing the 
system to find that the nonlinear fast wave equation contains a source term 
dependent on transverse gradients in the Alfven wave amplitude, hence show- 
ing that nonlinear mode conversion is permitted in their scenario. This is then 
demonstrated in a numerical simulation, where an initial condition perturbs the 
medium to generate an Alfven pulse with a profile oc sech^ (y) (where B = y). 
The pulse subsequently generates fast magnetoacoustic waves as the simulation 
evolves. They conclude that Alfven waves not only heat directly via the phase 



mixing dissipation method as per Heyvaerts & Priest ( 1983 ) but also indirectly 



through nonlinear coupling to the fast wave, which is itself dissipative. 
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Botha et al. (2000) extended this work and found that the efRciency of the 



mode conversion is determined by the frequency and amphtude of the Alfven 
wave, and by the gradient in the background Alfven speed. Their work showed 
that the fast wave component eventuaUy reached saturation, and concluded that 
(if the Alfven wave amplitude amplitude is high enough) the nonlinear mode 
conversion can be a significant sink of Alfven wave energy that otherwise would 
contribute to the linear phase mixing damping mechanisms as per Heyvaerts & 



Priest (1983). McLaughlin et al. (2011 ) further extended the work into the visco- 



resistive case, where they found that the equilibrium density profile (and hence 
the location of heating) is significantly modified by the the visco-resistivity and 
the ponderomotive effects (i.e. drifting of the heating layer can occur). 

In this paper, we focus on the possibility of the nonlinear Lorentz force as an 
agent for mode conversion in the general 2.5D MHD scenario, and also consider 
how ponderomotive effects identified in 1.5D models carry over. To do so, we first 



extend the source term analysis of Nakariakov et al. (1997) to a MHD scenario 



which may be inhomogeneous in both magnetic induction and density (i.e. the 
background Alfven speed is inhomogeneous). Then, we emulate the numerical 
experiment of Nakariakov et al. ( ,1997, ) in light of our extended analysis to discuss 
previously unreported ponderomotive effects. 

Thus, the paper is presented as follows; in f|2]we derive wave equations for 
a general 2.5D cold plasma, to determine the set of nonlinear interactions per- 
mitted between the MHD wave modes. In f|3j we analyse nonlinear source terms 
forced by general Alfven waves, and further consider the case for harmonic Alfven 
waves in |3.1| &: [3?2] . In f|4]we present results of the numerical experiment. Finally, 
we interpret our results and draw conclusions in fj5] 



2. MHD wave equations in a general, 2.5D f3 — medium 

We first determine the role nonlinear terms of the MHD equations play in 
facilitating mode coupling in a ideal, 2.5D, /3 = plasma permeated by a 
general, potential, magnetic field, which we take as Bq = [Bx,By,0] such that 
V X Bo = 0, with a background density po — po{x,y). We set d/dz = 
throughout, i.e. take z as the invariant direction. The governing, ideal, j3 — Q 
nonlinear MHD equations are: 



9v 

~dt 



+ (v • V)v 



V X B 



X B 



~dt 
dp 
dt 



V X (v X B) 
-V • (pv) 



(1) 



where the standard MHD notation applies. Initially, the medium is flow-free 
(vq = 0), and is perturbed by finite amplitude perturbations of the form B = 
Bo + h{x,y,z,t), p = po+ pi{x,y,z,t), v = + vi{x,y,z,t). This yields: 



^-^(Vxb)xBo = ^(Vxb)xb-^^-fl + ^ 
dt p.po ppa Po dt \ Pa 



(vi • V)vi 
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dh 



- V X (vi X Bo) = V X (vi X b) 



dp 
~dt 



r + V • (povi) = -V • (pivi) 



(2) 



where the left-hand-side/right-hand-side governs the linear /nonhnear behaviour 
respectively. 

We now decompose into a;?/2-components such that vi = {vx,Vy,Vz) and 
h = {bx,by^^z), giving: 



dvz 
dt 



+ 



B„ 



dby dbx 



dt upo \ dx dy 



dVy 



Bx ( dby dbx 



dt 
1 



/ipo V dx dy 



d d 
+ By — ] bz 
UPo \ dx dy^ 

dh 
dt 

-dt^dx - 



d_ 
dy 



{VxBy - VyBx) 



dbz^ 
dt 



B, 



d 



d_ 
^dy^ 

dpi d , , d , , 

+ 7^ {pOVx) + [PoVy) 



dx 



dt dx dy 
where the nonlinear components Ni,. . . , N-j are: 



iVi 

N4 
Nr 



(3) 



N2 = 

Ns = 

Ni = 

N5 = 



1 

Pi dv, 

Po dt 

1 



^ dbx 
dy 



b ^+bz^ 

^ dx ^ dx 



1 



Pi 

Po 



d 



d_ 
dy 



b ^ 
dx 



Pl dVy _ 

Po dt 
1 d 

b-r 



1 + 



, dbx ^ , dbz 

bx -7^ 1- ^ 

dy dy 

Pi\ f d , d 



Po 
d 



1+^ 
Po 



ppo \ dx 
Pl dVz _ 
Po dt 

^ (Vxby - Vybx) 

d 

-g^i^xby-Vybx) 



""^d-x^'^dy 



d d 
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d d 

Note that the nonlinear Lorentz force in the momentum equation has been 
considered in the ahernative tension-pressure form 

(V X b) X b = (b • V) b - (b • b) 

so that the contributions of the different aspects may be compared. Here, 

(b • V) b is the nonhnear magnetic tension and V (b • b) /2 is the nonhnear 

magnetic- pressure gradient. 

We now seek wave equations governing the evolution of the permitted wave 
modes. To do so, we differentiate the momentum equation terms to link the 
system of equations. We then find the velocity components corresponding to the 
fast, slow (absent in /3 = 0) and Alfven wave by using the coordinate system 



of Thurgood & McLaughhn (2012). They use this system to not only derive 
governing wave equations, but to allow the formation of initial conditions that 
correspond to pure, distinct wave modes and to decompose propagating waves 
into constituent modes (we utilise this approach in our simulations in fj3| . In our 
2.5D, Cartesian analysis the velocity components arc thus ti_L = Vi • (z x Bq) = 
—ByVx + BxVy (fast), Vz (Alfven), and = vi - Bo = B^Vx + ByVy (longitudinal 
perturbations, where the slow mode is absent in /3 = 0). Hence, the governing 
equations are: 



"92 




92 Y 






^9y2j_ 



9^ 
9i2 



1 



B 



9 



9 



^ dx ^ dy 



9iV3 
dt 



1 



B 



9 



B, 



9 



^ dx ^ dy 



(5) 



92^11 



B. 



B, 



dN2 
dt 



(6) 



where ca = \/ Bq/ fipo is the background/equilibrium Alfven speed that varies in 
the xy-plane, such that ca = CA{x,y). Note that by setting the nonlinear terms 
of equations (4][6) to be zero (i.e. the right- hand-side) we revert to the linear 
regime, and can see that the linear, /3 = fast waves and Alfven waves are 
completely decoupled, and there are no disturbances along the magnetic field. 
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3. Source Term analysis 



We now consider conditions on the fluid variables that correspond to a pure, 
linear Alfven wave (as per Alfven ,1942^ i.e. waves driven only by magnetic 
tension) at an early time when other modes are taken to be absent (i.e. coupling 
has not yet occurred) to determine which terms of equations (4]|6) will act as 
sources of other modes of oscillation in the general MHD case (i.e. terms that 
contribute to the acceleration of velocity components that are initially zero in 
the absence of the corresponding wave mode), i.e. we perform a source term 
analysis. 

To do so, we take 7^ and ^ with — Vy = ^ by = pi = (and 
so v± = U|| = 0), as perturbations in the z-direction correspond linearly to a 
pure Alfven wave (and the linear Alfven wave does not perturb mass density). 
Note that this sets A^^ = for i = 3, . . . , 7, and simplifies A^i and The wave 
equations (4][6| become: 



9V 

9*2 



1 d 
MPo dt 



By 7^ 1 br, 

ox oy ^ 



(7) 



B, 



dy 



(8) 



9*2 



1 d 
fipo dt 



B. 



d_ 
dx 



B„ 



dy 



b. 



(9) 



again, the right-hand-side contains the nonlinear terms. 

Here, we see that there is a source term associated with vj_, thus in this case 
it is possible that d^v±/dt'^ ^ 0. Inspection reveals that equation ([t]) can be 
rewritten as: 

2 



9*2 



1 d 
ppo dt 



z X Bo • V 



(10) 



Hence, the Alfven wave can cause a nonlinear magnetic-pressure gradient to 
arise which in turn causes the excitation of fast magnetoacoustic waves, where 
gradients transverse to Bq are responsible. Note that the term z x Bp • V is 
the gradient in the direction across fieldlines. This concurs with the analysis of 



Nakariakov et al. (19971, which considered nonlinear effects in a unidirectional 



field and found coupling terms associated with transverse gradients (note that 
if we specifically consider equations [4] and [5] for such a scenario, we recover their 
wave equations, see Appendix |Aj) Hence, the Alfven wave exerts a transverse, 
nonlinear magnetic-pressure gradient in any situation in which it assumes non- 
zero gradients across the magnetic field. However as highlighted by Verwichte 



( 1999 ) the net perturbation will only be non-zero if the ponderomotive force 



averaged over a wave-period is non-zero. 

Similarly, we find that the Alfven wave can cause a nonlinear pressure gradient 
to arise which causes longitudinal perturbations to the equilibrium field. By 
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rewriting ^ as 



d^vii 1 d 



dt^ fipo dt 



Bo- V 



(11) 



we see that this is due to longitudinal gradients in the Alfven wave amplitude, 
i.e. this is the longitudinal effect of nonlinear magnetic pressure. 



We note that equations (10) and (11) can be integrated with respect to time 
(see Appendix [C|) to yield 



dv± 
~dt 



Z X Bq ^ 


(f 


fJ-Po 






f) 


fJ-Po V 





(12) 



i.e., a familiar force equation ('F = ma') dependent on the square of the ampli- 
tude of the Alfven wave field, conforming to the general definition of a 'pondero- 
motive force' discussed in §1. Here, wc have adopted the notation V_l = zxBq-V 
and V|| = Bg • V - these terms are the gradients transverse and longitudinal 
relative to the equilibrium magnetic field in the a;y-plane. 

One can also perform a source term analysis for fluid variables corresponding 
to an initially pure fast wave. Doing so, we find that the ponderomotive effects 
of the propagating fast wave do not facilitate coupling to the Alfven mode, the 
fast wave does not interact with the Alfven wave on any level, linear or nonlinear 
in a medium with an invariant direction (see Appendix [b|). Hence, we further 
can conclude that ponderomotive conversion in /? = is a one-way process from 
the Alfven to the fast magnetoacoustic mode. 

3.1. Harmonic Alfven wave 

To further investigate the nature of ponderomotive coupling, we reconsider the 
source terms of the fast wave and longitudinal equations when forced by a 
harmonic linear Alfven wave of form 

V, ^ Acose ^^{Ae'^) , e^ujt-k-r , b^=Six,y)v, 

where S{x, y) is a spatial scaling term between velocity and magnetic field per- 
turbations that is associated with the background Alfven speed (for a linear 
Alfven wave, b = zLy/JIpov). If we consider that linear Alfven waves propagate 
along magnetic field lines only, then the direction of the wavevector must be 
k = Bq. Wc also can derive a dispersion relationship from equation (jsj linear 
terms only) that links frequency, wave number and Alfven speed: = w^/fc^. 
Thus, we can consider the wavevector to be k = fck = fcBo = uj'Bq/ca, hence 
the function 9 = uj (^t — Bq ■ y/ca^ ■ Inserting the harmonic form of such a wave, 

the source terms driving the fast and longitudinal modes (equations ([t] and [9]) 
become: 



9t2 fipQ c\ 



(Bo • r) {ca) - caVi_ (Bq • r) ] cos 26 (14) 
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2 

r [(Bo • r) V|| (ca) ~ c^V,, (Bq • r)] cos 20 (15) 



Hence, equations (14) and (15) govern the transverse and longitudinal pertur- 
bations forced via the action of a harmonic Alfven wave. 

The forced motion is driven at double the frequency (cos 26) and is of order 
0{A^), i.e. key features of the ponderomotive effect previously reported in the 



literature (see e.g., Nakariakov et al. 1997 Botha et al. 2000 and McLaughlin et 
al. 2011 ). We find that perturbations to v± and uy are dependent on two distinct 
classes of terms, one associated with gradients in the Alfven-speed profile and 
the other with gradients in the geometry of the propagating pulse relative to the 
equilibrium magnetic field. 

3.2. Instantaneous, geometric terms: Vj^ (Bq • r) & Vii (Bq ■ 



It is known in 1.5D homogeneous MHD, for an Alfven pulse propagating along 
a fieldline, that due to gradients in the wave intensity (i.e. ponderomotive ef- 
fect) the leading flank acts to longitudinally-accelerate the plasma and the rear 
flank acts to longitudinally-decelerate the plasma, thus sustaining a cospatial, 
instantaneous perturbation that is transported along the magnetic field which, 
depending on the specific pulse geometry, may or may not cause a non-zero (net) 
longitudinal perturbation as it passes through the medium (i.e. the 'ponderomo- 
tive wings' of Verwichte [1999J) . Regardless of whether a net perturbation to the 
plasma occurs, this longitudinal daughter disturbance is continually sustained 
by the propagating Alfven wave and always remains cospatial to its progenitor 
(viz. the daughter occupies the same spatial region as the Alfven wave). If we 
impose conditions equivalent to 1.5D homogeneous MHD (i.e. equations 14 and 
where Vj^ — and Vii {ca) = 0) only one driving term remains, 



15 



V|| Bo-r 



i.e. this term governs the longitudinal daughter disturbance. 

Our equations indicate that, if we extend our consideration to the homoge- 
neous 2.5D case (i.e. allow nonzero transverse gradients), a transverse equivalent 
to Verwichte's ponderomotive wings, a transverse daughter disturbance, can exist 
when there are transverse gradients in the pulse profile across the equilibrium 
magnetic field, i.e.: 



Bn 



Such a profile, in homogeneous MHD, is imposed as an initial condition and 
would be maintained throughout, however in inhomogneous MHD a profile with 
transverse gradients can developss naturally via phase mixing, i.e where the 
pulse enters a region with a transverse Alfven-speed profile, thus assuming a 
pulse geometry with transverse gradients. We stress that such a transverse 
inhomogeneity of the Alfven-speed profile (i.e. a phase mixing region) can be 
dependent on variable Bq and poj and is not just upon a density inhomogeneity 
in a unidirectional field. 
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Figure 1. The Alfvcn-speed profile C4 (red) and the density profile po (blue). There is 
steep, continuous, transverse gradient in Alfven speeds (a phase mixing region) in the region 
-1 <x <1. 



3.3. Inhomogeneity terms: Vj^ (ca) & V|| (c^) 

Now permitting V(ca) 7^ and considering 2.5D inhomogeneous MHD, our 
equations show the ponderomotive effects not only depend upon the instan- 
taneous magnetic pressure perturbations of assumed pulse geometry discussed 
previously, but also upon magnetic pressure perturbations as the pulse geometry 



is altered from instant to instant (i.e. the full set of equations 14 and 15 1. 
This term is thus the ponderomotive effect due to Alfven-speed profile inho- 
mogeneities, where its transverse manifestation (which can be thought of as a 
phase mixing term) excites fast waves, and the longitudinal manifestation excites 
longitudinal perturbations (which can be thought of as a longitudinal dispersion 
term). If the net longitudinal perturbation is non-zero, due to the absence of gas 
pressure gradients in /3 = the perturbation will be static. 

Thus, the analysis of |3.ip.3| shows that the harmonic Alfven wave nonlin- 
early interacts with the medium via the longitudinal and transverse manifesta- 
tions of two classes of ponderomotive effect, namely the geometric effect (cospa- 
tial, ponderomotive daughter disturbances) and the V(cyi) effect (ponderomotive 
acceleration due to inhomogeneous Alfven-speed profile). Both may yield non- 
zero net perturbations of the medium, causing coupling to the longitudinal 
and transverse (fast magnetoacoustic) modes. As the two terms may interfere 
(constructively or destructively), the precise dynamics of ponderomotive mode 
conversion will vary on a case by case basis. 



4. Numerical Demonstration 

We now demonstrate the ponderomotive effects identified in jp] |3.3| by consider- 
ing simulations of a simple scenario , namely that of an initially uniform pulse 
separating in a unidirectional field stratified by a smooth, yet steep, transverse 
density profile (i.e. a transverse Alfven-speed profile). Such a scenario was con- 
sidered in Nakariakov et al. (1997), we emulate their results to demonstrate 
features that they did not report upon. We solve the nondimensionalised, non- 



SOLA: Verl_Accepted_Thurgood_McLaughlin_SolPhys.t6x; 19 February 2013; 2:15; p. lOJ 



Ponderomotive effects induced by Alfven waves in 2.5D MHD 



linear MHD equations using LARE2D (see Arber et al. 2001 for details of code 



and, e.g., Thurgood & McLaughlin 2012 for details of the nondimensionalisation 



procedure, noting that all quantities are nondimensional in this section) for the 
equilibrium magnetic field Bo = y where the plasma is structured by a density 
inhomogeneity of the form: 

, , 0.2 + e-4^ 

in the numerical domain x, y G [—8, +8] with a resolution of 1920^ grid points 
for a cold plasma initially at rest (vq = 0). 

The resulting (nondimensionalised) Alfven-speed profile, shown in Figure [l] 
is thus ca ~ 1 for x < ~1, ca ^ \/5 for a; > +1, with a steep but continuous 
transition between the two limits in the region — 1 < x < 1. In this scenario 
(due to the unidirectional field) the longitudinal direction here is simply y 
and the transverse direction is x. Thus the transverse and longitudinal velocity 
components of the simulation (vx and Vy) relates to the transverse and lon- 
gitudinal velocity components in the previous analysis {v± and wy) such that 
vj_ = -ByVx + BxVy = -Vx and = B^v^ + ByVy = Vy. 

We perturb the system by imposing an initial condition of the form 



= 2Asech^ 



with pulse-width parameter a — 0.25, which creates an Alfven wave (with apex 
at ?; = 0) which separates into two oppositely-travelling pulses with amplitude 
A — 0.001. The initial amplitude is taken to be small so that if a fast wave 
is nonlinearly generated, it will not disturb the equilibrium field through shock 
wave formation. Simple reflecting boundary conditions are employed, thus the 
simulation is halted once the pulses reach the boundaries. Figure [2] shows the re- 
sulting separation and propagation of the Alfven waves (seen in 1;^). Throughout 
time, the pulse geometry maintains its initial longitudinal (oc sech^y) and trans- 
verse (= 0) profiles outside of the phase mixing region (x < — 1 and x > -1-1). In 
the region, due to phase mixing, the profiles change in time, with the creation of 
transverse gradients in pulse geometry which become increasingly steep as time 
progresses. 

We now consider perturbations to the longitudinal velocity component Vy, 
shown in FigurejSj Here, we note two nonlinear features. The first is a disturbance 
that appears to propagate as the Alfven waves (i.e. along the fieldlines at the 
Alfven speed) yet is observed in the longitudinal velocity component. These 
are not independently propagating waves (in our /? = simulation there is no 
gas pressure to facilitate longitudinal oscillations) but instantaneous perturba- 
tions/disturbances that are sustained and carried by the propagating Alfven 
waves in Vz, i.e. these are the longitudinal daughter disturbances identified in 
§3.2, and is associated with longitudinal gradients in the pulse geometry (for the 
harmonic Alfven wave, this was of the form V||Bo -r). The longitudinal daughter 
disturbance sustained by the upwardly propagating Alfven wave is of the same 
sign, whereas the daughter sustained by the downwardly propagating wave is the 
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t = 0.000 



t = 0.6?^ 



t = 1 .87 







Figure 2. Contour plots of v^, showing the propagation of the Alfven wave, which undergoes 
phase mixing due to the transverse inhomogeneity in Alfvcn-speed profile. 




t = j,l 2b 




Figure 3. The longitudinal velocity component Vy for the same times as Figure[2] We find a 
nonlinear static perturbation localised at the initial position of the Alfven pulse and nonlinear 
cospatial disturbances that are transported with the Alfven waves. The colour bar is scaled by 
a factor of xlO^, i.e. by xA'^. 
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t = 0.000 t = 0.62b t=i.2S0 




u b 



Figure 4. The transverse velocity component v^: for the same times as Figure|2] Propagating 
fast waves emanate from the phase mixing region. As the Alfven wave undergoes phase mix- 
ing, independently propagating fast waves are generated and a cospatial transverse daughter 
disturbance develops. The colour bar is scaled by a factor of xlO^, i.e. by xA^. 



opposite to its progenitor. Their maximum amplitudes are Vy = ±4.999 x 10 ^ 
respectively, which is of O (0.5A^). We note that amplitude of the longitudinal 
daughter is weaker by a factor of approximately \/b in the higher Alfven speed 
region compared to its counterpart in the lower Alfven speed region 

The second longitudinal ponderomotive effect is a static perturbation as- 
suming the geometry of the initial pulse, which is also is of smaller amplitude 
in the higher Alfven speed region. The fluid velocity is perturbed such that 
Vy = ±3.031 X lO"'^ at X = -8 and Vy = ±6.678 x 10"^ at a; = +8 (in the vicinity 
oi y = 0), hence this perturbation is also weaker in the higher Alfven speed 
region by a factor of approximately \/5, i-c. the amplitude of the longitudinal 
perturbations is inversely proportional to ca- The perturbation is accompanied 
by a stationary mass density enhancement (pi), initially generated at the same 
order as the velocity perturbations (the same magnitude and with the inverse 
proportionality to ca) then grows linearly in time during the simulation. These 
features are the 2.5D analogue of the cross-ponderomotive effect noted in the 
homogeneous 1.5D MHD study of Verwichte et al. ( p99| , a non-zero longitudinal 
plasma perturbation caused by crossing pulses (which corresponds to our initial 
condition) . 

Figure |4] shows the transverse velocity component Vx , which is the velocity 
component associated with the fast magnetoacoustic wave. As reported and 
thoroughly detailed in Nakariakov et al. (1997), we find nonlinear disturbances 
in Vx are generated in the phase mixing region and propagate outwards from 
regions of high-to-low Alfven speed, independently of the Alfven wave, clear 
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evidence that the Alfven wave nonhnearly generates independently propagating 
fast magnetoacoustic waves as it undergoes phasemixing. As the Alfven wave is 
small, the fast waves saturates at a low amplitude {vx = —3.338 x 10~^) relative 
to its progenitor, in agreement with Botha et al. (2000). Aside from outwardly 
propagating fast waves of Nakariakov et al. (1997), Figure |4] reveals that as the 
Alfven wave undergoes phase mixing, a transverse cospatial disturbance arises 
in the region, i.e. we can see the transverse daughter disturbance. This is most 
clearly visible when comparing the later panels of Figure [4] to the corresponding 
times on Figure [2] (to determine the region cospatial to the Alfven wave). 



5. Conclusion 

In this paper, we present two main results: 



The analysis and conclusions of Nakariakov et al. (1997) extend to the 
general 2.5D MHD case, such that the nonlinear magnetic pressure exerted 
by an Alfven wave propagating through a region of variable Alfven speed 
facilitates mode conversion from the Alfven to magnetoacoustic modes in 
many MHD scenarios (to the fast mode via phase mixing and to the longi- 
tudinal mode via dispersion along iieldlines). 

The ponderomotive effect not only generates independently propagating 
magnetoacoustic waves but also longitudinal and transverse daughter dis- 
turbances, which remain cospatial and dependent upon the progenitor wave. 



After deriving wave equations for a general MHD system (equations 4][6) 



with a non-rotational invariance (i.e d/dz = 0), the source term analysis of 
S|3] demonstrates that the nonlinear magnetic-pressure gradients generated by 
a propagating Alfven wave can facilitate mode conversion from the Alfven to 
magnetoacoustic modes, providing that, the pulse assumes a geometry which 
yields a non-zero net force over the wave's period. A similar source term analysis 
for a propagating fast wave is presented in Appendix [B] which confirms the 
additional result that the process is one-way (i.e. the nonlinear Lorentz force 
does not permit conversion from the fast magnetoacoustic to the Alfven mode). 
This, of course, is entirely intuitive as pressure gradients (magnetic or otherwise) 
cannot act in the Alfven wave direction, as it is impossible to have gradients in 
an invariant direction (by definition). 

In §3.1 we consider the terms forced by a harmonic Alfven wave, and find that 
perturbations to v± and wy are generated at double the frequency and at the 
square of the driving/initial amplitude, in agreement with results reported in the 
literature. Such a result is necessary for ponderomotive effects, as the nonlinear 
magnetic-pressure gradient acts upon the square of (and hence the square of 



Vz). The forced wave equations (14 1 and (15) are found to contain two distinct 



sets of terms, each with a longitudinal and transverse manifestation, that are: 

• Geometric effects 

(a) V|| ^Bo • - Longitudinal daughter disturbances 
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(b) Vj^ (^Bo • r j - Transverse daughter disturbances 

• V {ca) effects 

(c) Vj[ (ca) - Perturbations due to longitudinal dispersion 

(d) Vj^ (ca) - Perturbations due to phase mixing 

Reducing to homogeneous 1.5D MHD, we find that one term forcing lon- 
gitudinal perturbations remains, corresponding to V|| ^Bq ■ above, which 
is dependent on longitudinal gradients in a function of pulse geometry and 
background magnetic field. This term is responsible for the cospatial pondero- 



motive wing reported Verwichte et al. (1999). Upon permitting a transverse 
direction (i.e. homogeneous 2.5D MHD), we find that a transverse analogue is 
permitted. Hence, to distinguish the two we refer to the longitudinal/transverse 
manifestations as the longitudinal daughter disturbance and transverse daughter 
disturbance respectively. When we further consider the inhomogeneous scenario 
and permit gradients in the Alfven speed, a second class of forcing term is per- 
mitted, again with transverse and longitudinal manifestations, which accelerates 
the medium as the Alfven wave passes through regions of inhomogeneity. The 
transverse manifestiation occurs where the Alfven wave undergoes phase mixing 
(regions with transverse gradients in Alfven speed) and the longitudinal manifes- 
tation occurs where longitudinal dispersion occurs (in regions with longitudinal 
gradients in Alfven speed), where we stress that the Alfven-speed profiles vary 
with Bq and po (i.e., phase mixing can occur in regions of uniform density, if the 
magnetic field has transverse stratification, an example of this is magnetic null 



point configurations, see e.g., Fruit & Craig 2006) 



To illustrate the phenomena implied by our analysis and interpretation of 



31-^3.3 we considered numerical simulations of plasma with a unidirectional 



field structured by transverse density profile (the same scenario considered in 



Nakariakov et al. 1997). In addition to the results detailed in Nakariakov et al. 
(1997) - namely, the generation of independently propagating fast waves, we 
find clear evidence for the existence of the longitudinal daughter disturbance, 
which remains cospatial to its progenitor throughout its transit. We note that 
the disturbance is differs by a factor of a/S in the higher Alfven speed region 
(relative to its amplitude in the lower Alfven speed region). An explanation is 



found by reconsidering equation (15) in the limit x — — oo and x — +oo for 
our scenario. Specifically, for the scenario there is no longitudinal dispersion 
(V||C^ = 0) and both the frequency and the longitudinal geometric gradient is 

identical at both extremes (thus we arbitrarily take a; = 1 and V|| ^Bq • = 1 

without loss of generality), and the scaling function is taken as — fipQ. The 
equation implies 

1 

2c^a' 

Given that in the simulation ca {x = — oo) = 1, ca {x = +oo) = 5 with the 
Alfven wave amplitude A — 0.001 we find that iiii {x = — oo) = 5.000 x 10^^ 



max « - — 
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and I'll (a; = +00) — 2.236 x 10^^ which is in excellent agreement with the 
reported amplitudes of the longitudinal daughter disturbances in the homoge- 
neous regions. Hence, the longitudinal daughter disturbances of Alfven waves 
are generated with amplitude of O (A^/2c^). 

The simulations also reveal the transverse manifestation of the cospatial 
disturbance, the transverse daughter. This is observed to develop over time (ini- 
tially, transverse gradients in pulse geometry zero) as the Alfven wave undergoes 
phase mixing, remaining cospatial to the phase mixed part of the Alfven wave. 
To our knowledge, this is the first time this phenomena has been reported. 

An additional effect noted in the simulation is that of the static perturbation 
aty — that is manifest in the longitudinal, y-components of fluid variables (e.g. 
in Vy, Figure [3|. We believe that this is the much the same cross-ponderomotive 
effect as reported in Verwichte et al. (1999), despite the fact that our scenario is 
inhomogeneous 2.5D as opposed to the homogeneous 1.5D analysis they perform. 
This is because (a) the cross-ponderomotive effect will not have any transverse 
action as Alfven waves cannot separate in the transverse direction (and there is 
nothing to suggest such a phenomena in the simulation data) , and (b) there is no 
inhomogeneity along the field lines, and hence dispersion terms (VyCyi) will not, 
in the specific case, impact upon the separation dynamics. As our simulation is 
inhomogeneous (with transverse stratification), it highlights that perturbations 
due to the cross-ponderomotive effect are (like the other ponderomotive effects 
observed) inversely proportional to the speed ca- It is however likely that in 
general, dispersion during pulse crossings will yield a more complicated effect. 
For cases where propagating Alfven wave pulses are likely to meet (in particular 
where reflection may occur) and the scenario is sufficiently nonlinear it would be 
necessary to investigate the ponderomotive effects on interaction between Alfven 
waves further. We suggest that it would be necessary integrate equation ( 15 ) 
over a separation period of a general wave form on Vz (from which the form of 
bz follows) that corresponds to the general solution satisfying the variable-speed 
1.5D wave equation uu — c{y)uyy for travelling pulses from a Gaussian-like initial 
condition (i.e. a variable speed D'Alembert solution). Such an analytic solution 
can be found by transforming the variable speed wave equation to a constant 
coefficient Klein-Gordon equation, however the profile of the wave speed has to 
conform to various conditions (of which, our scenario in §4 does not satisfy). 
See Grimshaw et al. ( j2010j ) for a comprehensive overview of the homogenisation 
process. 

During the separation of the pulse, we also note the generation of a stationary 
density enhancement (initially of O (A^/2ca)) that subsequently grows linearly 
in time. Again, this is in agreement with that reported by Verwichte et al. ( 1999 ). 
Linear density instabilities were shown to be a general feature oi P — MHD 
by Falle and Hartquist (2002). There, the authors analyse the eigenmodes of 
cold MHD and show that there is a Jordan Mode Instability, associated with the 
absent /zero-speed slow waves, which causes such density enhancements. This 
instability is to shown arise when the parallel (i.e. slow) velocity component is 
not constant. The cross-ponderomotive effect of a separating Alfven wave (or two 
crossing Alfven waves), as observed in our simulations, excites such a mode by 
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acellerating the parallel velocity component by applying a (nonlinear) magnetic 
pressure gradient. 

Throughout the paper we have considered the cold plasma regime. If the /3 ^ 
MHD wave equations (equations 4][6 if gas pressure gradients are permitted) are 
considered under the conditions for an initially pure Alfven wave as per ^ we 
find that gas pressure gradients do not contribute, i.e. the /3 ^ equivalents 
of the source term derived equations ([7| - ( 15 ) are unchanged, i.e. the way in 
which a passing Alfven wave nonlinearly perturbs the medium is identical to 
that detailed in this paper. Hence, where non-zero longitudinal perturbations of 
the medium occur (which are static in /? = 0, e.g. the cross-ponderomotive 
perturbation in Figure [s]), gas pressure gradients will subsequently arise to 
transport the disturbance longitudinally, i.e. slow magnetoacoustic waves will be 
generated. Thus, in /? = 0, ponderomotive mode conversion is permitted between 
the Alfven wave and both magnetoacoustic modes. In addition to the implications 
for mode conversion in an extension to a /3 7^ regime, we note that the cross- 
ponderomotive effect, and subsequent excitation of a Jordan mode instability, 
will still generate large density enhancements in low-/? plasmas as per Falle and 
Hartquist ( 2002 1 . Thus, crossing Alfven waves will directly contribute to the 
creation of inhomogeneity in low-/3 plasmas, due to their cross-ponderomotive 
effects. 

We note that aspects of mode conversion, coupling and interaction, explored 
in this paper from the perspective of applied nonlinear forces, can be investi- 
gated using alternative methods, in particular a MHD instability approach. The 
relationship between the results presented here, and literature on Alfven wave 



instabilities (such as, e.g., Derby 1978 Goldstein 1978 Champeaux et al. 1997 
Webb et al. [2001 ) will be explored in future work. 

We conclude the paper by highlighting that all of the analysis can be re- 
peated with different types of invariant coordinate system. For instance, we could 
consider azimuthal invariance, or the general coordinate system of Thurgood 
& McLaughlin (20121 for zero-helicity topologies, and yield the same results. 
Given that an invariant coordinate is a necessary condition for the existence of 
true Alfven waves as per Alfven (1942) (see, e.g., Parker ,1991, and section on 
coordinate systems in Thurgood & McLaughlin 2012 their section 2.3.1), we 
believe our results and conclusions extend to any scenario in which an Alfven 
wave (a wave driven by magnetic tension only) can exist. We have shown that a 
nonlinear Lorentz force of a propagating Alfven wave can generate independently 
propagating fast and slow magnetoacoustic waves, dependent on the interplay 
between Alfen speed profile and pulse geometry (viz., the form of the daughter 
disturbances), which varies on a case by case basis. As the transient properties 
of the magnetoacoustic modes are fundamentally different to the Alfven wave, 
such conversion facilitates the indirect transport and dissipation of (initially) 
Alfven wave energy to plasma regions that are inaccessible in the linear MHD 
regime. In modelling wave behaviour in sufficiently nonlinear solar plasmas, the 
ponderomotive effects of propagating waves must be evaluated. 
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Appendix 



A. Recovery of wave equations for a medium with a homogeneous 
magnetic field and a transverse density profile (Nakariakov et al. 
1997) 

Here we demonstrate that our general wave equations for 2.5D MHD (|4] & [s]) 
reduce to those considered in the case of Nakariakov et al. (1997), which con- 



sidered a unidirectional homogeneous magnetic field structured by a transverse 
density profile. To do so, we take the field as constant in y, i.e Bq = Bqy, 
Bx — and By = Bq and consider density as a function transverse to tlie field, 
i.e. Pa = po{x). Thus, the Alfven and fast wave equations (|4]&[5]) become 



c\{x) 



Bo dNe 



dt fJ-Poix) dy 



(16) 



r a2 



4 



9^ 



-Bo 



dNi 



dy 



c\{x) 



cm 

dx 



(17) 



The Alfven wave equation ( 16 ) is only superficially different from the equation 



governing the Alfven wave in Nakariakov et al. ( 19971, which (in their notation) 
is 



c\{x) 



1 f dN2 
Po{x) \ dt 



Bo dN5 
4tt dz 



(18) 



where the apparent difference is due to 



The analysis of Nakariakov et al. (1997) considered y as the invariant 



direction, hence Vy corresponds to the Alfven wave (as opposed to Vz in 
this paper as z is invariant). Additionally, d/dy corresponds to d/dz in 
this paper and Nakariakov et al. ( 1997 ) respectively. 

Hence, the nonlinear terms have a different ordering; their TVi, N^, N4, 
N^, Nq ,Nj corresponds to our Ni, N3, N2, N4^, Nq, N^, respectively. 
The nonlinear terms originating in the equation of motion (iVi, iV2, N3) are 
not exactly identical. They differ by a factor of Pq^ (we divided through 
the equation of motion by po to group it with the nonlinear terms, they left 
it with the linear terms, e.g., our iVi is equal to their Ni/po). 



Nakariakov et al. (1997) use cgs units, as opposed to SI in this paper. 



The fast wave equation of Nakariakov et al. (|1997l is 

Bo dNi 



dt^ 



d^ 
dy'^ 



1 



dNi 



Bo dNe 



A-K dz A-K dx 



(19) 



the above only differs from our fast wave equation as per the aforementioned 
superficial differences. Additionally, the nonlinear terms (right hand side) differ 
by a factor of —Bo- The explanation is simple- in this case uj_ = Vi • z x Bq — 
—BoVx - Both correspond to a transverse perturbation of the field in this scenario 
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and hence wave equations on either describe the evolution of the Alfven wave; 
in a homogeneous field it is simpler to use v^, for an inhomogeneous field v± 
is required. Expressing v± in terms of Vj. and simplifying reduces our fast wave 



equation (17) to the comparable form 



dNi ^ Bo dN^ Bo dN^ 
dt fJ.poix) dy fJ-Poix) dx 



(20) 



B. Nonlinear effects of a propagating fast wave. 



We can set = and bz = and place restrictions on v^, Vy, b^, and by 
that correspond to ~ &|| ~ 0, v± 0, b± and pi 7^ to consider the 
ponderomotive effects caused by a fast wave, yielding: 



5y2 



^dN^ ON, fdN, 
Ox "Jt; Oy h Cj^ 



dt 



1^ 



= 



dt 



V dx 



(21) 
(22) 



1 d 
fipo dt 

1 d 
PPo dt 



Bxby 



dbx 
dy 



dby 
dx 



Bybx 



dby 
dx 



dbx 
dy 



[dbx 


dbyV 


\dy 


dx J 



(23) 



where subscripts T and P indicate whether a term is contributed by magnetic 
tension or magnetic pressure respectively, and b± — h ■ z, x Bq is the transverse 



perturbation of the magnetic field. The key result here is that equation (22) 
shows that that the fast wave does not interact with the Alfven wave on any 
level; linear or nonlinear. 



However, equation (23 1 contains source terms, thus a propagating fast wave 
does cause nonlinear field-aligned disturbances that will cause, in certain cir- 
cumstances, nonlinear Lorentz force coupling (via both nonlinear magnetic pres- 
sure/ponderomotive force and nonlinear magnetic tension ). Departing from the 
/? = case, these source terms remain, i.e. the magnetoacoustic modes are 
coupled linearly by gas-pressure gradients, and nonlinearly by the Lorentz force. 



Equation (21 1 shows that such an initial fast wave undergoes self interaction. 
In this case, this is due to a combination of nonlinear Lorentz force, convective 
acceleration (contained within iVi and N2 which are permitted to be non-zero) 
and nonlinear induction {N4 and Nq). Interestingly, the nonlinear induction term 
can be rewritten as 



4 



aiV. dN4 



dx 



dy 



-4 



d^ 

dx"^ 



dy'^ 



V • z X b 



i.e. the nonlinear induction term is effectively of the same form as that of the 
linear motion due to velocity perturbations across the equilibrium field, however 
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is instead motion due to perturbations across the induced magnetic field (though 
at the equihbrium, not induced, Alfven speed). 

FinaUy, we note that although there arc no slow waves in the /3 = case, 
the nonlinear source term due to the Lorentz force acting upon the Alfven wave 
(right hand side equation 22 ) would remain unchanged for /3 7^ 0. This source 
term remains zero if an Alfven wave is initially absent. Thus, the nonlinear 
Lorentz force does not facilitate slow to Alfven mode conversion either. Hence, 
the analysis shows that for any 2.5D MHD scenario that that there is a one-way, 
nonlinear mode conversion mechanism between the Alfven and (both) magne- 
toacoustic modes, that will become manifest when a phase mixed/ dispersive 
Alfven wave (via propagating in the vicinity of a inhomogeneous Alfven-speed 
profile) assumes a geometry that exerts a non-zero average force over its period 
perturbing along and across the magnetic field. 



C. Derivation of force equations (12) and (13) 



We demonstrate that equations (10 1 and (11 



time to yield the force equations ( 12 ) and 



I can be integrated with respect to 
131 (i.e., with zero value 'integra- 



tion constants') by simply considering an alternate, direct derivation from the 
momentum equation. 

Under the source-term conditions of fjs] {vz ^ and 6^ 7^ with Vx ^ Vy = 
bx ~ by = pi = 0) the X and y components of the momentum equation simplify 
to yield 

^=_J_b — = 

dt upo ^ dx 2/xpo dx 



dv. 



y_ 



I 



db. 



MPo ^ dy 



1 dbl 
2ppa dy 



From which equations (12) and (13) readily follow by constructing equations in 
terms of vj_ and v\\ from vj^ = vi • (z x Bq) = —ByVx + BxVy and v\\ = vi • Bq = 



ByVy. 
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